A model description of the supersolid state in YBCO. 
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I employ a semiphenomenological model introduced in A. M. Tsvelik, A. V. Chubukov, Phys. Rev. 
Lett. 98, 237001 (2007) to describe the state with co-existing superconductivity (SC) and charge 
density wave (CDW) recently discovered in YBCO. The SC and the CDW order parameter fields are 
united in a single pseudospin and can be rotated into each other. It is suggested that disorder creates 
isolated pseudospins which become centers of inelastic scattering of electrons. It is suggested that 
this scattering is responsible for the logarithmic upturn in the resistivity p{T) ~ — InT observed at 
low doping. 
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Recent x-ray measurements in clean YB2Cu306+a; 
show that in a certain doping interval superconduc- 
tivity in this material coexists either with Charge Den- 
sity Wave (CDW) order or at least with strong CDW 
fluctuations. The CDW peak emerges at temperatures 
below 150K and is not accompanied by enhancement of 
the magnetic fluctuations. According to 3], an incom- 
mensurate SDW appears at smaller dopings, closer to the 
ones where the superconductivity disappears altogether. 
Application of magnetic field suppresses the supercon- 
ductivity, but the CDW remains. Its presence should 
lead to a reconstruction of the Fermi surface (FS) which 
has been invoked by Sebastian et.al. [1] as an explana- 
tion of the results of de Haas van Alphen effect mea- 
surements [3, [5], [6]. However, since the bulk of the FS 
is already gapped at these temperatures, as is evident 
from the sharp drop in the integrated intensity of spin 
excitations occuring at ~ 20 meV Q, the CDW cannot 
be driven by such quasiparticle mechanisms as the FS 
nesting. It is reasonable therefore to suggest that low 
energy physics of underdoped cuprates is dominated by 
collective excitations and the quasiparticles play a sec- 
ondary role in determining their phase diagram (not in 
the transport where their role is important). 

This is the point of view adopted in this paper. I con- 
sider a model which describes coexistence and competi- 
tion of SC (superconductivity) and CDW in a state of 
preformed pairs where quasiparticles play a secondary 
role. At low temperatures the pairs may either local- 
ize (the CDW insulator) or exist in a supersolid state 
where charge density fluctuations coexist with supercon- 
ductivity. The bulk of the bare FS is supposed to be 
gapped; there is a small quasiparticle FS whose area 
scales with doping. Within this model at zero temper- 
ature CDW-l-SC(supersolid) phaseQ is separated from 
the pure CDW by a Quantum Critical Point (QCP). It 
would be logical to associate with it the line of QCP's in 
H — X plane {H is magnetic field, x is doping) observed 
in the conductivity measurements where SC has 

been suppressed by magnetic field. However, the present 
x-ray data [1] contradict this suggestion: the CDW peak 
disappears below =0.09, well before the critical doping 



=0.04. 



Universality of Quantum Critical Point physics relieves 
one from concerns about microscopic mechanisms be- 
hind the adopted low energy effective theory giving a 
certain freedom in its choice. I choose a model capable 
to describe a competition and coexistence between CDW 
and SC and broadly compatible with the known cuprate 
physics. 



The phenomenological model I use to extract the low 
energy properties are inspired by ID physics, but can be 
justified by the DMFT results. These results as described 
in fTl]| suggest that, as far as observables are concerned, 
strongly interacting 2D systems and ID systems are quite 
similar. First and foremost, the DMFT suggests that the 
pseudogap phenomenon is not necessarily accompanied 
by any kind of order [l^ which is similar to dynami- 
cal generation of spectral gaps in ID models. Second, 
there is a striking similarity in detail between results of 
[lH,[l3l and the physics of doped fermionic ladders as 
discussed, for instance, in fl3|. Such similarity has in- 
spired now famous Yang- Rice- Zhang (YRZ) ansatz [l^ 
successfully employed as a phenomenological tool in the 
cuprate physics. In the latter case there is also a spin gap 
and singularities in the electronic self energy not accom- 
panied by any obvious density wave order. The obvious 
difference though is that in ID systems all these phenom- 
ena occur at any bare couplings due to the always present 
nesting, but in 2D one needs strong interactions. 



Although I am mostly interested in universal low en- 
ergy physics, it would be desirable to use as a starting 
point more or less realistic lattice models. To illustrate 
the physics relevant to the QCP I choose a model of 
crossed one-dimensional chains forming a grid introduced 
in (l6j . The grid network is supposed to emulate the in- 
terplay between CDW (or stripes) and superconductivity 
in terms of preformed Cooper pairs, which appear nat- 
urally in ID if the interaction is attractive. Each chain 
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contains electrons with a strong local attraction: 
i?. - {-t. [C+ (r + x)a (r) + H.c 

r=(aonx,) 

- UN^{r)N^{r)-f,J2Nair)}, 



(1) 



where N — J^a ^t^a- The same Hamiltonian describes 
chains running parallel to y-axis. I denote its creation 
and annihilation operators B^,B. The electrons are al- 
lowed to tunnel from the chains along x to the chains 
along y and viceversa at the intersections. The Hamilto- 
nian of crossed chains is 

H = YH.,[ny\+Y^ Hy K] + (2) 

71 y 7lx 

t± ^ [C^iaon^r. aony)B^{aoNn:r,0'oNny) + H.c], 



where oq is a distance between parallel chains. I assume 
that tj_ « tx « U. 




FIG. 1: The grid-like lattice of model ((2}. The hollow cir- 
cles mark lattice sites where B, and C, C*" operators are 
defined. The solid and dashed lines mark tunneling pathes 
along y and x axes respectively. 

As I said, the sole purpose of model ([2]) is to surve as a 
seed for universal low energy physics which will be sensi- 
tive only to its gross features so that all nonrealistic ones 
will be "forgotten". To obtain this low energy descrip- 
tion I will employ a sequence of stepwise transformations. 
A validity of each subsequent transformation will be con- 
fined to smaller and smaller energy. The first step on this 
path is to transform the Hubbard model ([T]) to a model 
of magnet. For energies well below the single particle gap 
~ U model ^ is equivalent to the isotropic spin S=l/2 
Heisenberg model in magnetic field which role is played 
by chemical potential /i: 
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- ^(7'^(i" + x)T-(r) +H.C.] + 



^(r + xK(r)] +2^iT^{v)], 



(3) 



where r° are Pauli matrix operators, = f^/U. T^(r) 
create and destroy electron pairs located at site r and 



stands for the deviation local charge density from 1. This 
transformation valid for the attractive Hubbard model is 
similar to a more familiar transformation for the repulsive 
case where low energy degrees of freedom are described 
by real spin operators. One important difference is that 
the exchange integral of the transverse spin components 
in the attractive case is ferromagnetic; this signifies that 
the pairs carry zero momentum. The exchange interac- 
tion in model (jSj can be made anisotropic by including 
longer range charge interactions. The repulsive ones gen- 
erate an Ising-like anisotropy and as a consequence will 
open a spectral gap at = (a commensurate CDW), 
but this gap is removed ones the consentration of pairs 
deviates from 1/2. The coupling between the perpendic- 
ular chains has the same Heisenberg form, but with a 
smaller exchange integral J± — t\/U . 

The next step in the derivation of the low energy action 
proceeds by bosonization of model ([3]): 
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T^{x = nax) — kp/n + ^ —d^Q 
sm{2kFX + %/2^e)(-l)" + .. 



(4) 



(5) 



where dots stand for less relevant operators and kp — 
TTTiQ. Here the mutually dual bosonic phase fields $ and 
Q are governed by Hamiltonian ^ eplained below. The 
bosonized form of the Hamiltonian is valid as yet lower 
energies that the range of validity of ([3]) , namely, below 
Jx- I remind the reader that no stands for concentration 
of pairs. In the real YBCO where the preformed pairs 
coexist with the quasiparticles a relationship between kp 
and the concentration of carriers is not straightforward. 
The magnitude of kp may be even doping independent 
as suggested by the experiments 0- 

The presence of the constant term in ^ implies that 
at nonzero doping each chain experiences a Zeeman- 
like field with wave vectors 2Tm/Nao. Therefore, when 
■K — 2kpa{) = 21:11/ N , one of the harmonics of the peri- 
odic potential strongly couples to the staggered magneti- 
zation of the T-field (charge density wave field for actual 
fermions) ([5]). The continuum limit Hamiltonian valid at 
energies below ^ tl./U is 
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Asin(%/27re„ + Qa;)|, 



(7) 



where ['^{x'),dxQ{x)] = i5{x — x'), g± ^ J±,X ^ kp 
and Q = ^(1 - 2n/N) - 2kp. K ^ 1 for the large U 
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Hubbard model, but can be renormalized by longer range 
interactions. The velocity is v ^ J. 

At the commensurate point Q = the spectrum of 
each individual chain is gapped (assuming K < A). The 
field on a given chain is pinned by the periodic poten- 
tial generated by the chains crossing the given one. This 
is a CDW ordered insulating state where the interchain 
coupling J± is irrelevant. Therefore in order to have an 
effect this coupling must exceed some critical value. 

Let us concentrate on the interchain coupling. The 
operator 

* = ei^* (8) 

taking part in the interchain coupling describes a charged 
bosonic field of charge 2e (the Cooperon). As for any 
Lorentz invariant theory the dispersion of particles in 
the sine-Gordon model ([7]) has a relativistic form: e = 
■\/ (vq)^ + . For K ~ \ the spectrum consists of 
a degenerate triplet with m ~ M and a singlet with 
m — \/3M, where M ^ X^/i'^-^) is determined by pa- 



rameters of model ([T]). Two components of the triplet 
(soliton and anti-soliton) are charged Q = ±2e. Thus 
Cooper pairs in the insulating state are gapped excita- 
tions which gap is directly related to the CDW ampli- 
tude. They are also coherent ones since operator (|8]) 
has a non-zero matrix element between the vacuum and 
a single-soliton state [l3|. The corresponding Green's 
function is (see jl8i] ) 

= ico + io/-^vqr-M^ + 

where the dots stand for the incoherent part, M is the 
soliton mass and Z is a numerical factor which exact 
value is of no importance here. Th coherent part of ^ 
can be represented as a correlation function of complex 
relativistic massive Bose field (Cooperon). Taking into 
account the interchain interaction in RPA I obtain the 
following Cooperon Green's function for coupled crossed 
chains: 
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where g±_ = Zg±. The excitation spectrum of the 
Cooperons is 



u;' = M2 + gV2 ± ^(Afgl)2 + cos2(2a), 

tana ^ qy/q^. (11) 

Close to critical point g±'^ — M, at v\q\ << M, we have 

uj^_ ^ M{M - g~_i_) + v^q^/2 + 0{q^). (12) 

Hence at energies << g~± the action can be approximated 
by the Lorentz invariant $^ theory in (2-|-l)-dimensions: 

where r = M{M - g\) and u ~ ]VP. At g~j_ < M 
the system is an insulator, the conductivity of pairs is 
a ~ exp(— Y^/T), at g~± > M it is a supersolid. It is 
remarkable that at the critical point the sheet resistivity 
of pairs is 19]: 

R{lj) = (/i/l-037(4e2))p(cj/T) (14) 

with p(0) = 1 and h/{2ef « 6A5kn. 

By construction of model ^ the CDW order always in- 
cludes two density waves with {Qo,0) and (0,Qo) (Qo = 
ir/ao — 2kp) wave vectors. Since these CDW's are driven 



by the periodic potentials of the crossed chains, they are 
not a result of phase transition. With increase of temper- 
ature their magnitude simply gradually diminishes from 
its T = value (sin(v^e^,j,)) - M^/^ ~ A^/t^-f^). 
Since in the superconducting phase the spectral gap is 
closed only in at a specific point in momentum space, 
namely at q = 0, the CDW amplitude, which magnitude 
is related to the average gap, in that phase does not van- 
ish. The characteristic temperatures at which the CDW 
order is noticable are of the order of M. On the other 
hand, the superconducting order occurs at temperatures 
determined by the stiffness and hence related to the mean 
field value of l^fj: Tc ~ —rv^/u and hence sufficiently 
close to the QCP, where M » Tc, the CDW order pre- 
dates the superconductivity. This conclusion agrees with 
the observations in YBCO [l|. It is also important that, 
due to the breaking of translational invariance, the su- 
perconducting stiffness is not directly related to the den- 
sity of pairs, but is non-universal being determined by 
parameters of model p^ . 

To make contact with real cuprates two additional fac- 
tors should be considered; quasiparticles and disorder. 
The presence of quasiparticles in the normal state of un- 
derdoped YBCO is well attested by the ARPES measure- 
ments done in zero magnetic field |2l| and by the de Haas 
van Alphen measurements done in strong magnetic field. 
The existing difficulties in figuring out an exact shape 
and evolution of their Fermi surface are not relevant for 
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the present discussion. Whatever their FS is, the quasi- 
particles can be easily accomodated in the present model, 
as it was discussed in [l^. The action of the preformed 
pairs (jl3|) should be augmented by the action of quasi- 
particles. As I said in the beginning the thermodynam- 
ics is dominated by the collective modes, but the normal 
state transport is dominated by the quasiparticles. This 
is due the fact that in the normal state the conductiv- 
ity of the Cooperons cannot exceed certain value which 
order of magnitude is determined by (|14l) . The latter 
one is the minimal metallic conductivity and in a clean 
YBCO th quasiparticles short circuit the conductivity of 
the Cooperons. 

In the region of phase diagram where CDW was ob- 
served, YBCO is a very clean compound as shown by the 
mobility measurements Q and NMR ^2^]. The inverse 
life time of these quasiparticles in the strong magnetic 
field regime is estimated as ~ 5meV [l^. A disorder 
introduced in CuO planes causes two effects: the Tc is 
suppressed and p(T) reveals InT upturn at low temper- 
atures. This strong logarithmic upturn has not been ex- 
plained by standard localization theories. According to 
[iCf] , the coefficient at In T is proportional to the impurity 
concentration suggesting that it originates from a strong 
scattering on isolated centers. This brings to mind Kondo 
effect as a possible explanation (as it was suggested in 
p^). However, the conventional Kondo effect where the 
inelastic scattering originates from localized spins is in- 
consistent with the fact that the In T upturn seems to be 
impervious to strong magnetic fields applied to suppress 
the superconductivity. I suggest instead a charge Kondo 
effect where the role of local spins is played by singlet 
local pairs created by disorder [2^. It is essential that 
the disorder is not sufficiently strong to localize quasipar- 
ticles (otherwise there would be no Kondo effect). Exis- 
tence of delocalized quasiparticles in the relevant range 
of disorder and magnetic fields is guaranteed by the fact 
that the In T upturn has been observed in samples whose 
resistance slightly above the transition was smaller than 
the quasiparticle localization shreshold: R{Tc) < 25.4fcl7. 

Let us now imagine that we are in the CDW regime 
where in a clean sample pair excitations would have a 
finite energy. However, the disorder by radically weaken- 
ing the interaction of a given pseudospin with its neigh- 
bours can create a localized two-level state with zero en- 
ergy which strongly scatters the quasiparticles. Such lo- 
calized states are described as pseudospins as in model 
([3]). An isolated pseudospin is free to flip in the pro- 
cess of interactions with the quasiparticles (not shown in 
(HI), as I said above to describe the transport one has to 
add quasiparticles) which causes their inelastic scatter- 



ing. The Kondo Hamiltonian describing the interaction 
between pseudospins and quasiparticles is given by 



Hk - ^6(fc)a+(fc)a,(fc) + 2^M(r)r3(r) 

R 

y~i ^ ei('^-'^'>{[7±(k-k')T+(r)a^(k)a;(k') + i?.. 

k,k',r 

+7||r3(r)a+(k)a,(k')}, (15) 

where a J , Qo- are creation and annihilation operators of 
quasiparticles. The c?-wave symmetry of the order param- 
eter wave function is accomodated in the ^-dependence of 
the coupling constants. A strong contribution to electron 
scattering comes from the pseudospins whose chemical 
potential is smaller than the Kondo temperature Tr- (its 
value depends on the coupling constants and the quasi- 
particle bandwidth). Such pseudospins will dominate the 
resistivity provided the distribution function of fi(r) has 
a sharp peak at small /i's. 

In this paper I have shown that the theory formulated 
in brings together such features of low energy physics 
of the underdoped YBCO as coexistence of superconduc- 
tivity and charge density wave and the InT upturn in 
the resistivity. The latter effect is caused by disorder and 
quasiparticles. It is quite possible that the model has a 
more general applicability since there are indications of 
low doping QCP in other cuprates. 

The model postulates existence of intermediate energy 
scale g± (smaller than the pseudogap) serving as the up- 
per cut-off for the theory In [16|] it was associated 
with the emergence of enhanced diamagnetism; it may 
also be related to the emergence of CDW observed by x- 
ray scattering As I have already mentioned within 
the present setting this order emerges gradually and is 
not accompanied by any thermodynamic anomalies. The 
change of sign of the Hall coefficient jE],[2^] occuring in 
the the same temperature range may also originate from 
the quasiparticle Fermi surface reconstruction caused by 
the CDW. ARPES experiments on BiSCO also point to- 
wards existence of such intermediate energy scale below 
the pseudogap psj . 
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